ETATS DE COMPACITE MAXIMALE POUR LES MELANGES 
BINAIRES DE GRAINS SPHÉRIQUES : ÉTUDE PAR SIMULATION 

NUMÉRIQUE 



Jean- Noël Roux*"), François Chevoir( a) et Fabrice Toussaint^ 

(a) Laboratoire des Matériaux et des Structures du Génie Civifl Institut Navier 
2 Allée Kepler, 77420 Champs-sur-Marne 

(b) Lafarge, Centre de Recherches 
95 rue du Montmurier BP 15 
38291 Saint-Quentin Fallavier 

RÉSUMÉ : Nous appliquons la simulation numérique (par la méthode des éléments discrets, également 
appelée dynamique moléculaire) à l'étude des assemblages désordonnés de compacité maximale de par- 
ticules sphériques, monodispersées ou bien en mélanges binaires, le rapport des diamètres étant fixé à 
trois. La simulation conduit à une définition pratique de la compacité maximale des systèmes désor- 
donnés, comme celle des configurations mécaniques d'équilibre sous chargement isotrope de particules 
rigides non frottantes. Une étude statistique des effets de taille d'échantillons permet d'évaluer, dans la 
limite des grands systèmes, les compacités maximales des assemblages monodispersés, pour lesquels on 
retrouve la valeur classique 0.639, et bidispersés, pour deux valeurs distinctes de la fraction volumique de 
grosses particules. Les calculs mettent en évidence l'influence sur l'état final et sa compacité de la durée 
de l'étape initiale de l'assemblage, dans laquelle les grains sont agités et interagissent par collisions. L'al- 
longement de cette phase d'agitation entraîne une évolution graduelle de l'état du système dans laquelle 
l'augmentation de la compacité s'accompagne de la nucléation d'un ordre cristallin dans le cas des billes 
de même taille et d'une ségrégation progressive selon la taille dans le cas des mélanges binaires. Même 
si cette évolution est lente et son effet probablement peu sensible dans nombre de situations pratiques, 
il convient de définir la compacité maximale des états désordonnés comme celle que l'on obtient dans 
la limite des assemblages rapides. Nous indiquons pour finir quelques prolongements possibles de cette 
étude préliminaire de la compacité des mélanges. 

1 Introduction 

De nombreux matériaux du génie civil sont composés d'un assemblage de grains solides, 
seuls (sols pulvérulents, graves non traitées, granulats, ballasts), en présence d'un fluide (sols 
saturés, ciments frais, matériaux de chaussées lors de leur mise en place), ou inclus dans une 
matrice solide (bétons, enrobés bitumineux en service). La compacité - c'est-à-dire la fraction $ 
du volume de matériau occupée par les particules solides - est toujours au nombre des grandeurs 
utilisées pour caractériser l'état de ces matériaux, pour tenter de prédire leur comportement ou 
aider à leur formulation. Souvent on se réfère à une compacité maximale $*, implicitement sup- 
posée liée à la seule géométrie des grains. Ainsi, pour la formulation des ciments 0}, on corrèle 
la viscosité d'une suspension à la compacité maximale des grains solides. En mécanique des 
sols, où on utilise souvent l'indice des vides e = au lieu de la compacité, des procédures 



'Laboratoire mixte, rattaché au LCPC, à l'ENPC et au CNRS. 



1 



normalisées en donnent des valeurs maximales et minimales, tandis que l'étude du compor- 
tement mécanique des matériaux granulaires souligne l'intérêt de la notion de compacité (ou 
d'indice des vides) critique, qui est approchée dans les états d'écoulement plastique en grandes 
déformations 0. 

Les méthodes de simulation discrète [3], qui connaissent actuellement un grand essor, donnent 
accès à l'ensemble des degrés de liberté d'un assemblage de grains solides. Le présent article a 
pour but de montrer qu'elles permettent de préciser la notion d'état de compacité maximale, en 
dégageant les influences de la géométrie et de la micromécanique, et de déterminer $* pour des 
systèmes simples dont la granulométrie n'est pas trop étalée. 

Pour un rapport des diamètres établi à trois, on a choisi d'étudier des mélanges avec deux 
proportions différentes p de volume solide correspondant aux grosses particules : p = 0.5 et 
p = 0.7, et de limiter le nombre total N de particules à quelques milliers dans les échantillons 
simulés. 

Le bilan de cette étude préliminaire, avec les résultats obtenus, l'influence de N et les fluc- 
tuations statistiques entre échantillons différents, est présenté et brièvement discuté dans le §131 
Auparavant, le §[2 donne quelques définitions utiles et présente rapidement les méthodes numé- 
riques employées, tandis que des résultats connus dans le cas du système monodispersé sont 
rappelées au §[21 Nous terminons en suggérant (^3) quelques perspectives quant à la poursuite 
de l'application de la simulation numérique discrète aux assemblages granulaires polydispersés. 

2 Définitions, méthodes numériques. 

2.1 Qu'est-ce qu'un état de compacité maximale ? 

Avec des billes de la même taille, la compacité maximale de l'assemblage aléatoire est 
classiquement estimée à une valeur légèrement inférieure à 0.64 [@][5]|. 

La notion d'état désordonné de compacité maximale est a priori un peu floue. Il s'agit en 
principe d'un maximum géométrique, concernant des particules rigides et impénétrables, mais 
qui exclut certaines configurations considérées comme trop ordonnées (en particulier celles dans 
lesquelles on peut clairement déceler un arrangement cristallin, au moins localement). 

Une remarque permet de préciser cette définition : si l'on se donne une définition de la com- 
pacité $ adaptée aux conditions aux limites choisies, les configurations qui réalisent un maxi- 
mum local de $, sous les contraintes d'impénétrabilité, coïncident exactement avec les confi- 
gurations d'équilibre stable de grains rigides non frottants soumis à une pression isotrope 0. 
«Maximum local » signifie qu'il n'existe pas de mouvement augmentant $ qui ne viole pas 
les contraintes stériques. Le problème posé en termes géométriques a donc naturellement un 
sens mécanique : partout où les contraintes d'impénétrabilité sont actives, apparaît une force 
normale entre les deux grains en contact, qui est un multiplicateur de Lagrange associé à cette 
contrainte. 

Il est donc tout naturel de chercher à obtenir ces états en soumettant des grains non frottants 
à une pression isotrope. C'est un procédé fréquemment employé en simulation numérique pour 
obtenir des configurations denses [r7l[8ll9l[TÔl. 

Un état de compacité maximale est un état d'équilibre d'en ensemble de grains rigides sans 
frottement dans les contacts, soumis à une pression isotrope. 
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2.2 Qu'est-ce qu'un état désordonné de compacité maximale ? 

On dispose de différents indicateurs (paramètres d'ordre) pour la cristallisation des systèmes 
de « sphères dures » ifTTl [T2l ITÔ1 . Celle-ci se manifeste à partir d'une compacité d'environ 
0.49 en présence d'une agitation thermique maintenue (dans le modèle de fluide pour lequel 
les collisions conservent l'énergie cinétique), mais elle est relativement facile à éviter dans 
la pratique des calculs si on se place à $ > 0.55 [fTT|. Avec des billes macroscopiques qui 
dissipent l'énergie dans leurs contacts et échappent à l'agitation thermique spontanée, on peut 
retrouver la tendance à la cristallisation lorsque l'assemblage est soumis à des vibrations ou à 
un cisaillement oscillant |[T3l . 

Le phénomène de cristallisation disparaît lorsque la largeur de la distribution des diamètres 
(écart-type d'une distribution gaussienne rapporté à la valeur moyenne) est de l'ordre de 6% du 
diamètre moyen |[T4||. 

Ce qu'il faut craindre alors, c'est la ségrégation des particules par taille, autre cause possible 
d'inhomogénéité. Les mélanges de sphères de différentes tailles, si on les considère comme des 
molécules d'un fluide ou des particules colloïdales soumises au mouvement brownien, peuvent 
exhiber un diagramme de phase assez complexe, qui dépend fortement du rapport des tailles, 
comporte éventuellement plusieurs phases fluides et plusieurs phases solides, etc.. On pourra 
par exemple consulter la réf. lfî~5l . et la copieuse littérature à laquelle renvoie cette étude. A 
densité élevée, la structure la plus stable implique toujours une séparation de deux phases de 
compositions différentes. Comme dans le cas monodispersé, en revenant aux particules macro- 
scopiques (c'est-à-dire sans agitation thermique ou mouvement brownien), on observe que la 
manipulation d'un mélange de particules de tailles différentes (qu'il s'agisse d'une vibration 
continue, de secousses répétées séparées par des temps de repos, etc..) tend à induire le même 
type d'organisation que l'agitation thermique spontanée des systèmes moléculaires ou colloï- 
daux : on assiste à une ségrégation progressive tfT6l[T7l[T8l . De même la mise en écoulement des 
matériaux granulaires conduit à la séparation des populations de grains selon leurs tailles |fî~9ll . 

Avec les sphères d'une seule taille on peut définir un état de compacité maximale aléatoire 
si on a un maximum local de $ sous les contraintes d'impénétrabilité dans l'espace des confi- 
gurations et si on ne décèle aucun germe de cristallisation locale. 

Il est logique d'adopter une définition analogue pour les systèmes bidispersés, en remplaçant 
la cristallisation par la ségrégation. 

On constate que, pour des échantillons assez grands, les états désordonnés de compacité 
maximale possèdent tous {approximativement, et non pas exactement, en raison des évolutions 
lentes vers la cristallisation ou la ségrégation) la même compacité $*, qui est celle que nous 
étudions ici. 

Les résultats rappelés ci-dessous dans le cas des billes d'une seule taille et présentés au 
^4] pour les mélanges binaires indiqueront dans quelle mesure, pour le procédé d'assemblage 
indiqué au §12.31 ci-dessous, la compacité peut dépendre de la durée de la phase initiale dans 
laquelle les grains sont agités. 

Un état désordonné de compacité maximale peut être défini comme un état d'équilibre d'en 
ensemble de grains rigides sans frottement dans les contacts, soumis à une pression isotrope, 
pour lequel les évolutions vers l'ordre cristallin ou vers la ségrégation des mélanges restent 
négligeables. 
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2.3 Procédure de simulation 



Pour produire des états homogènes, commodes à caractériser et dépourvus d'effets de bord, 
on utilise des conditions aux limites périodiques. 

La méthode de calcul utilisée est la dynamique moléculaire, appellation abrégée en DM 
ci-dessous, (que l'on désigne souvent comme « méthode aux éléments discrets» ) avec des 
sphères assez rigides et des ingrédients d'inertie et de dissipation visqueuse appropriés - ces 
derniers garantissant une restitution nulle dans les collisions binaires. La raideur des contacts est 
telle que dans l'état final 1' «interpénétration » typique des sphères en contact (c'est-à-dire, bien 
sûr, la déflexion normale du contact) soit de l'ordre de quelques 10~ 5 <i, d étant leur diamètre. 
On trouvera des explications et des références plus détaillées sur la méthode dans un article 
récent 0. 

Au préalable, les particules sont disposées soit « au hasard » (un procédé qui impose de 
commencer les calculs dans une configuration de compacité assez faible), soit sur les sites 
d'un réseau régulier. Dans une seconde étape, la configuration de départ est mélangée : on 
donne aux grains une certaine énergie cinétique et on les laisse interagir par des collisions qui 
conservent cette énergie (chocs élastiques), à la manière des molécules d'un fluide (qui seraient 
toutefois dépourvues d'interaction à distance). Lorsque la position de départ était sur un réseau, 
on s'assure alors que cet ordre initial a bien disparu dans l'étape de mélange. 

La durée de la phase de mélange et son influence possible se mesurent au nombre de colli- 
sions moyen par grain, à une compacité donnée. 

Une possibilité intéressante, en particulier dans le cas de systèmes polydispersés, pour les- 
quels il peut être difficile de déterminer une configuration initiale au hasard avec une compacité 
qui ne soit pas trop faible, est l'utilisation d'un algorithme qui permette simultanément de mé- 
langer et de compacter, pour que le recours à la DM puisse n'intervenir que dans une étape 
plus tardive. C'est ce que fait l'algorithme de Lubachevsky et Stillinger 11201 I2T1 (abrégé en 
LS ci-dessous), qui consiste à calculer des collisions entre des objets qui se dilatent à vitesse 
constante, de sorte que la compacité augmente au cours du temps. 

Les résultats présentés au §|4]font tous appel à cette méthode dans une première étape, et ses 
effets sont testés ci-dessous dans le cas monodispersé. 

3 Résultats pour les billes monodispersées 
3.1 Procédure directe, temps caractéristiques 

On a pu vérifier que l'on obtenait avec un seul diamètre de billes la valeur $* = 0.637 ± 
0.001 en effectuant l'étape de mélange à $ = 0.45, (il suffit alors d'une dizaine de collisions 
par grain pour faire disparaître toute mémoire de la configuration initiale) puis en achevant le 
compactage par DM avec des particules déformables et des chocs mous. Cette valeur résulte de 
moyennes sur une série de 5 échantillons de N=4000 sphères. Une autre série de configurations 
de 1372 sphères donne 0.636 ± 0.0010. L'erreur due à la légère déformabilité des contacts était 
avec le modèle retenu estimée à A$ = 1.2 x 10 -4 (donc négligeable, vu l'incertitude sur $). 

2 On convient ici, et dans toute la suite, d' indiquer comme intervalle d' incertitude un écart-type de part et d' autre 
de la moyenne 
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Pour donner une idée approximative de la durée du processus de compaction (deuxième 

étape), qui est calculé à pression contrôlée P, on peut utiliser le temps caractéristique r = 

associé à l'inertie des grains de masse m et de diamètre d soumis à des forces d'ordre Pd 2 . La 
compacité finale est approchée (à moins de 0.01 près) en quelques dizaines de r. Ensuite, si 
on désire une configuration bien équilibrée (dans laquelle le réseau des contacts et les valeurs 
des forces sont déterminées avec une bonne précision), il faut poursuivre les calculs pendant 
un temps qui peut atteindre quelques centaines de r (il faut noter qu'il est relativement diffi- 
cile de stabiliser des grains non frottants dans une position d'équilibre). Pour avoir un point de 
comparaison avec une expérience de laboratoire, on peut considérer par exemple un processus 
de dépôt sous gravité. Une telle procédure est couramment employée dans les laboratoires de 
mécanique des sols, et en contrôlant débit et hauteur de chute libre - on parle alors de pluviation 
contrôlée ||22ll23l - elle constitue un procédé de fabrication reproductible d'échantillons homo- 
gènes dont la compacité peut être variée à la demande. Par la simulation numérique Il24ll25l . on 
peut comparer les assemblages de billes assemblées par pluviation aux configurations obtenues 
par les processus de compactage isotrope idéaux, dépourvus d'effets de bord, tels que celui que 
nous employons dans la présente étude. On vérifie alors que les états que nous étudions ici sont 
bien similaires à ceux que l'on peut trouver dans des situations plus proches d'une expérience 
réelle. 

Dans le procédé de dépôt gravitaire, la pression typique, celle qui règne dans la couche 
proche de la surface libre où les grains se réarrangent, correspond au poids d'une dizaine de 

couches granulaires. En prenant donc P ~ 10—^-, on aura r ~ 06 ^ onne envrron 

3 x 10~ 3 s avec des grains millimétriques. Les temps trouvés (de l'ordre du dizième de seconde 
pour l'assemblage par déposition d'une couche de 10 grains d'épaisseur, ou de la seconde pour 
remplir un récipient de 10cm de hauteur) semblent donc comparables entre le calcul numérique 
et la fabrication d'échantillons en laboratoire. 

Cette compacité pour les assemblages de sphères d'une seule taille est en bon accord avec 
les résultats d'expériences en laboratoire flU |5), et avec d'autres travaux numériques. Ainsi 
O'Hern et al. Il26ll par exemple, donnent $* = 0.639 ± 0.001, résultat obtenu par extrapolation 
vers la limite N — > oo. Ces auteurs ont mené à bien une étude systématique de l'influence de 
la taille du système sur la compacité maximale. Ils ont observé une légère influence du nombre 
N de sphères sur la valeur moyenne dans les échantillons de N sphères, ($)(iV) Celle-ci tend 
à augmenter vers sa valeur asymptotique $* dans le limite N — > oo comme : 

($)(iV) = $* - 5 A^" a , (1) 

avec ô = 0.12 ± 0.03 et a = 0.47 ± 0.05. Ils ont également déterminé l'influence de TV sur la 
largeur w (largeur à mi-hauteur d'une courbe en cloche) de la distribution des compacités dans 
les échantillons de taille N : 

w = w N~ Q , (2) 

avec w = 0.16 ± 0.04 et fi = 0.55 ± 0.03. Les équations Q] et [2] s'appliquent aux résultats 
de B2l pour N > 64. 



5 




3.2 Effet d'une agitation prolongée : algorithme LS 

Dans un second temps nous avons testé l'usage de l'algorithme LS pour ne recourir à la DM 
qu'à une densité plus élevée. 

Si l'on utilise l'algorithme LS jusqu'à obtenir $ = 0.632 (avec 6000 collisions par par- 
ticule en moyenne pour arriver à cette compacité) on constate que l'achèvement ultérieur du 
compactage par DM produit des compacités plus élevées que la valeur trouvée précédemment : 
$ = 0.6423 ± 0.0004. Il est difficile de repérer un ordre local dans ces nouvelles configurations, 
l'évolution vers l'arrangement cristallin restant très embryonnaire (ceci peut se contrôler quan- 
titativement avec des « paramètres d'ordre » comme dans lfi"Tl[T2lO . Mais si on prolonge encore 
l'étape de calcul par la méthode LS (comme cela a été fait dans Il27ll28l ). on peut obtenir toutes 
les compacités intermédiaires entre celles que nous avons citées et la valeur ~ 0.74 du 
réseau cristallin parfait, correspondant à des systèmes partiellement ordonnés selon les motifs 
cubique à faces centrées ou hexagonal compact. 

On voit donc qu'il n'est pas possible en principe de considérer la valeur de la compacité 
maximale $* des systèmes désordonnés comme unique. Toutefois, si l'usage de l'algorithme 
LS est poussé jusqu'à $ = 0.59 avec seulement 100 collisions par bille, on obtient à nouveau 
la valeur précédente ($ ~ 0.637 pour N = 4000, légèrement inférieure pour N = 1372). 

On retiendra donc que la compacité maximale la plus usuelle correspond à un assemblage 
rapide, et qu'on la retrouve en limitant la durée des phases d'agitation. 

Cette interprétation est confortée par l'absence ou la brièveté d'une telle phase de vibration 
dans les procédés pratiques d'assemblages tels que le dépôt sous gravité (pluviation) ou le 
damage par couches, ainsi que par les résultats déjà cités de 11261 . qui sont en très bon accord 
avec les valeurs de <3>* que nous avons obtenues en limitant la durée de l'agitation. En effet, 
la méthode numérique employée par O'Hern et al. 11261 est censée correspondre à la limite des 
compressions très rapides (l'énergie potentielle élastique, associée à la déformation des contacts 
intergranulaires dans une configuration à volume imposé, est directement minimisée par un 
algorithme de descente du gradient, et la limite où le volume choisi donne des déformations 
nulles fournit la valeur de la compacité cherchée). 

3.3 Rôle du frottement 

La procédure d'assemblage précédente, qui consiste à faire décroître le volume de la cellule 
(périodique) contenant les grains jusqu'à ce que le réseau des forces de contacts vienne équi- 
librer chaun des grains et contrebalancer la pression extérieurement appliquée, mise en œuvre 
en présence de frottement intergranulaire, conduit à des configurations mécaniquement équi- 
librées de compacité moindre IZllSlIâl- De plus, en présence de frottement intergranulaire, on 
doit également noter que les grains peuvent être assemblés, pour les mêmes contraintes, dans 
des états assez différents, ce qui est bien sûr conforme à l'expérience, puisque l'on caractérise 
classiquement l'état d'un échantillon de sable par sa compacité. 

La définition des états de compacité maximale que nous avons donnée plus haut, comme 
états d'équilibre sous chargement isotrope de grains non frottants, conduit donc à considérer les 
procédés de compactage (vibration, cisaillement, lubrification) comme destinés à limiter ou à 
circonvenir les effets du frottement dans les contacts. 

Pour des billes sphériques identiques, la compacité minimale que l'on peut atteindre dans 
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une configuration d'équilibre dépend de la valeur du coefficient de frottement, et divers procédés 
peuvent conduire à des valeurs différentes, les plus basses étant de l'ordre de 0,56 à 0,58 avec 
des coefficients de frottement compris entre 0,2 et 0,5. 

Signalons enfin que la simulation numérique a permis de constater que la compacité et 
le nombre de coordination d'un assemblage de billes frottantes pouvaient varier indépendam- 
ment BUEES! . 

3.4 Effet d'un cisaillement 

Jusqu'ici nous avons considéré exclusivement des configurations isotropes. Or, il est pos- 
sible d'appliquer à des assemblages de sphères sans frottement une contrainte de cisaillement, 
la configuration restant en équilibre mécanique, avec des forces de contact normales qui équi- 
librent les forces extérieurement appliquées. Nous avons ainsi observé que l'on pouvait appli- 
quer des contraintes principales a 1 > o 2 = cr 3 jusqu'à atteindre des rapports 01/03 de 1.2 à 
1.25, avec un modèle de sphères quasi-rigides. Lorsque le déviateur des contraintes est aug- 
menté, par paliers, à partir d'une configuration assemblée de façon isotrope, on constate que le 
système se déforme, car le réseau des contacts se réajuste pour pouvoir supporter les nouvelles 
contraintes. A l'équilibre on obtient une distribution des orientations des contacts anisotrope, 
avec un plus grand nombre de vecteurs normaux aux contacts parallèles à la direction principale 
majeure des contraintes. Il est intéressant de constater que les réarrangements se produisent sans 
changement notable de la compacité. Cette remarque nous permet donc de généraliser la défini- 
tion que nous avons adoptée pour un état de compacité maximale, en abandonnant la condition 
de contraintes isotropes : c'est simplement un état d'équilibre de grains rigides non frottants. 

Dans le langage de la mécanique des matériaux granulaires (ou des sols O), les assemblages 
de sphères non frottantes possèdent un angle de frottement interne (faible, entre 5 et 6 degrés, 
mais non nul), mais sont dénués de dilatance. Pour un assemblage de grains parfaitement lu- 
brifiés, et donc non frottant, la compacité $* que l'on évalue en imposant une pression isotrope 
semble coïncider avec la « compacité critique », celle d'un échantillon continuellement cisaillé 
en régime quasi-statique. Cette valeur $* pourrait donc en principe avoir une certaine pertinence 
pour l'étude de la rhéologie des pâtes, pour lesquelles on s'intéresse à la compacité associée à 
la divergence de la viscosité, en régime newtonien. Une étude expérimentale récente ÏÏ291 par 
une équipe du LMSGC montre néanmoins que la viscosité diverge pour des compacités infé- 
rieures à $*, et il semble que le comportement rhéologique des pâtes granulaires soit influencé 
par l'état de surface des particules solides, que l'on ne peut pas considérer comme parfaitement 
lubrifiées par le fluide. 

Par ailleurs, nous ignorons si les assemblages de sphères polydispersées non frottantes sont 
également dépourvus de dilatance. 

Enfin, il convient de rappeler que la compacité critique des grains secs dépend du coefficient 
de frottement et est inférieure à la compacité maximale (elle est d'environ 0,595 pour les billes 
identiques lorsque le coefficient de frottement de contact vaut 0,3 0). 
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4 Résultats pour les mélanges binaires. 



Nous avons simulé la compaction isotrope d'échantillons de billes de diamètre d et D = 3d. 
Le tableau Q] donne la liste des différentes séries d'échantillons ainsi préparées (il rappelle aussi 
les résultats obtenus au laboratoire avec les billes d'une seule taille). La figure Q] représente la 



Code 


P 


N 


N 2 


M 


$0 


*1 


Ucoll 


($*) 




A00 





4000 





5 


0.45 


0.45 


50 


0.637 


0.001 


A01 





4000 





5 


0.05 


0.632 


5000 


0.6423 


0.0004 


COO 





1372 





10 


0.45 


0.45 


50 


0.636 


0.0012 


COI 





1372 





8 


0.05 


0.59 


80 


0.636 


0.001 


B10 


0.5 


2916 


104 


5 


0.05 


0.68 


860 


0.7041 


3 x 10" 4 


Bll 


0.5 


2916 


104 


5 


0.05 


0.70 


4400 


0.7070 


1.1 x 10" 4 


CIO 


0.5 


1372 


49 


10 


0.05 


0.68 


850 


0.7034 


6.1 x 10" 4 


Cil 


0.5 


1372 


49 


10 


0.05 


0.67 


2200 


0.7061 


7 x 10~ 4 


D10 


0.5 


700 


25 


20 


0.05 


0.68 


850 


0.7034 


10~ 3 


B20 


0.7 


2916 


232 


5 


0.05 


0.69 


120 


0.7164 


10~ 3 


B21 


0.7 


2916 


232 


3 


0.05 


0.69 


1650 


0.7207 


4 x 10~ 4 


C20 


0.7 


1372 


109 


10 


0.05 


0.70 


180 


0.7165 


9.6 x 10~ 4 


D20 


0.7 


704 


56 


20 


0.05 


0.69 


270 


0.7165 


1.2 x 10~ a 


D21 


0.7 


704 


56 


20 


0.05 


0.69 


850 


0.7185 


1.05 x 10~ 3 


E20 


0.7 


352 


28 


40 


0.05 


0.69 


100 


0.7143 


2 x 10" 3 



Tab. 1 - Différentes séries d'échantillons. N est le nombre total de billes, p la proportion en 
volume de grosses billes, N 2 leur nombre, M le nombre d'échantillons. L'étape LS donne lieu 
à ricoii collisions par particule pour compacter de $ à $i, et la compacité finale obtenue, après 
l'étape de calcul par DM (à partir de $i) est ($*) en moyenne, l'écart-type étant A$*. 

vue en coupe d'un échantillon de 2916 billes (dont 232 grosses) dans un état de maximum de 
la compacité, dont la figure [2] donne une vue en perspective. Les résultats de compacité sont 
présentés sur la figure [3l en fonction de N~ l l 2 (qui donne des variations à peu près linéaires 
de la moyenne et de l'écart-type de $* selon les relations Q] et [2]). Les tailles d'échantillon 
étudiées sont tout aussi représentatives avec le mélange binaire que dans le cas monodispersé : 
on ne constate ni une dépendance en N, ni un niveau de fluctuations de $* sensiblement plus 
élevé dans le mélange. Suite à un réglage progressif des paramètres de l'étape L-S, le nombre 
de collisions par particule lors de cette phase d'agitation et gonflement est variable entre les 
différentes séries d'échantillons (tableauQ]). Toutefois, une tendance se dégage nettement : plus 
l'étape d'agitation se prolonge dans le processus d'assemblage, plus la compacité finale sera 
élevée - résultat encore analogue au cas monodispersé. 

Au-delà de la seule valeur de la compacité, la simulation numérique - en particulier lorsque 
l'on poursuit les calculs, comme nous l'avons fait ici, jusqu'à obtenir de bons états d'équilibre 
mécanique - fournit bien d'autres informations. On a ainsi accès au réseau des contacts et aux 
forces qu'ils transmettent. Le tableau [2] donne les différents nombres de coordination z^-, définis 
comme le nombre moyen de contacts d'une particule d'espèce i (i = 1 pour les petites, i = 2 
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FlG. 1 - Vue en coupe d'un échantillon avec N = 2916 et N 2 = 232. Les particules figurées 
en gris appartiennent à la cellule de base dans laquelle sont conduits les calculs, et qui se répète 
périodiquement. 

pour les grosses) avec des voisines d'espèce JE Une certaine proportion x des Ni = N — N 2 
petites billes ne transmet aucune force, alors que toutes les grosses sont sollicitées (c'est-à-dire 
que les contraintes d'exclusion stérique sont actives). D'autre part, le tableau [3]donne les valeurs 
en unités de Pd 2 , P étant la pression appliquée, des moyennes des forces F^ et des rapports 
Fjp entre écarts-types et moyennes (caractéristiques de la largeur de la distribution, ^y) pour 
les contacts entre une bille d'espèce i et une autre d'espèce j. Ces données concernant le réseau 
des contacts, conformément aux résultats pour la compacité, indiquent que les assemblages 
granulaires obtenus, hormis une légère influence de la durée de la phase d'agitation, peuvent 

3 On a bien sûr 

£21 = f£\ 3 

z 12 p \dj 
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FlG. 2 - Vue en perspective d'un échantillon avec N = 2916 et N 2 = 232. 



bien être considérés comme préparés dans le même état, dont l'échantillonnage statistique est 
satisfaisant, indépendamment de N dans l'intervalle étudié. Dans la limite des faibles pressions 
ou des billes rigides, on doit avoir [0 un nombre de coordination moyen z* , évalué avec les 
seules sphères qui portent des forces, égal à 6 : 

Ni N 2 

-(z n + z 12 ) + — wr(z 2 i + z 22 ) = 6, (3) 



N-x Ni ' N-xoN-t 

correspondant à un réseau de contacts isostatique. Les résultats du tableau [2] sont proches de 
satisfaire cette condition, ce qui montre que la géométrie d'un assemblage de grains très faible- 
ment déformables en contact est décrite avec précision. 

Le tableau |2] montre aussi que les échantillons préparés avec une plus longue étape d'agita- 
tion (ceux dont le code se termine par « 1 », voir le tableau [[]), et qui possèdent une compacité 
plus élevée, présentent apparemment certaines prémices de ségrégation : le nombre de coordi- 
nation z 22 est plus grand, signe que les grosses particules se rassemblent. Il conviendrait bien 
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0.73 



0.72 



0.71 



0.7 



~i 1 1 r 



~i 1 — j AUI 1 1 1 1 1 1 1 1 1 r 

p = : * + 0.09 (+ lent) 

p = : $ + 0.09 (rapide) 

J A0 



p = 0.7 



p = 0.5 



j i i i i i i i_ 



-I ï 

O'He 



BOO 



?r n et 



ai. $ 



+ O.09 



B21 



B20 



plus lent 

j C20 



Bll 



BIO 



rapide 



plus lent 

T Cil 
J CIO 

rapide 

j i j i L 







0.01 



0.02 
N-i/2 



0.03 



D21 
D2C 



J I L 



0.04 



FlG. 3 - Compacité maximale $* représentée fonction de 1/yN pour les valeurs p = 0.5 
et p — 0.7 étudiées (les barres d'erreurs s'étendent ici sur un écart-type de part et d'autre 
de la moyenne). Les codes indiqués sont ceux du tableau [H et on a distingué par la couleur 
bleue les calculs pour lesquels une phase d'agitation prolongée a augmenté la valeur de la 
compacité finale, par comparaison aux processus d'assemblages rapides. Pour faire figurer sur 
le même graphe les résultats sur les systèmes monodispersés, on a décalé les valeurs de $* 
correspondantes de 0.09 vers le haut. Les résultats de ll26ll . résumés par les relations © et ©, 
sont figurés en vert. 



sûr, au-delà de cette étude préliminaire, d'étudier ce phénomène de plus près, et d'examiner si 
la tendance peut être poussée plus avant. Ces mêmes échantillons se distinguent également par 
un rôle accru des contacts entre gros grains dans la transmission des forces : la valeur de F 2 2 
(tableau [3]) montre que ces contacts sont non seulement plus nombreux, mais aussi davantage 
sollicités. 

La force normale de contact fki entre le grain k et son voisin /, dans la présente approche, 
sont des grandeurs géométriques, des multiplicateurs de Lagrange associés à la contrainte d'im- 
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Code 


«il 




221 


^22 


z* 


x 


A00 










6.074 


1.3% 


B10 


4.11 ±0.03 


0.877 ±0.003 


23.7 ±0.8 


2.34 ±0.05 


6.068 


5.6% 


Bll 


4.08 ±0.04 


0.876 ±0.014 


23.7 ±0.1 


2.48 ±0.1 


6.064 


6.1% 


CIO 


4.08 ±0.06 


0.890 ±0.018 


24.0 ±0.5 


2.24 ±0.19 


6.061 


5.7% 


Cil 


4.12 ±0.04 


0.885 ±0.023 


23.9 ±0.6 


2.47 ±0.12 


6.089 


5.4% 


DIO 


4.09 ±0.05 


0.884 ±0.026 


23.9 ±0.7 


2.30 ±0.3 


6.058 


5.5% 


B20 


2.76 ±0.05 


1.35 ±0.015 


15.6 ±0.2 


3.33 ±0.05 


6.055 


13.6% 


B21 


2.65 ±0.05 


1.33 ±0.02 


15.3 ±0.2 


3.73 ±0.13 


6.053 


15.7% 


C20 


2.74 ±0.04 


1.37 ±0.02 


15.9 ±0.2 


3.44 ±0.10 


6.079 


13.6% 


D20 


2.62 ±0.07 


1.38 ±0.04 


15.9 ±0.5 


3.48 ±0.25 


6.045 


14.7% 


D21 


2.66 ±0.11 


1.37 ±0.04 


15.9 ±0.4 


3.53 ±0.23 


6.049 


14.4% 


E20 


2.70 ±0.15 


1.36 ±0.06 


16 ±0.7 


3.37 ±0.32 


6.043 


14.0% 



Tab. 2 - Nombres de coordinations, proportion a;o de petites billes inactives. Notations définies 
dans le texte. 



Code 


F u 


^12 


F22 


-Ml 


r 12 


r 22 


B10 


0.827 


1.24 


4.47 


1.32 


1.36 


1.14 


Bll 


0.780 


1.24 


4.88 


1.36 


1.41 


1.15 


C10 


0.823 


1.25 


4.51 


1.33 


1.38 


1.17 


Cil 


0.796 


1.21 


4.80 


1.34 


1.40 


1.15 


D10 


0.829 


1.24 


4.51 


1.33 


1.36 


1.15 


B20 


0.960 


1.43 


4.93 


1.45 


1.51 


1.18 


B21 


0.800 


1.30 


5.70 


1.55 


1.62 


1.19 


C20 


0.951 


1.41 


4.81 


1.44 


1.49 


1.19 


D20 


0.931 


1.38 


5.08 


1.45 


1.52 


1.17 


D21 


0.901 


1.36 


5.18 


1.45 


1.52 


1.18 


E20 


0.953 


1.39 


5.04 


1.44 


1.49 


1.18 



Tab. 3 - Forces moyennes, écart-type réduit des distributions de force, par familles de contacts. 
Notations définies dans le texte. 

pénétrabilité. Si le rayon du grain k subit une petite variation 5R, alors le volume du système 
tout entier, si la liste des contacts reste inchangée (ce qui est toujours possible dans un échan- 
tillon fini avec un SR assez petit), varie au premier ordre de 

5V = Y J 1^5R. (4) 

1 

De plus, le rapport fki/P dépend dans la limite des grains rigides de la seule géométrie du 
paquet de billes en contact, du fait de l'isostaticité de cette structure H. On conçoit donc que 
ce type de résultats puisse fournir des indications pour la recherche de granulométries optimales, 
l'équation © pouvant donner une idée de l'effet d'une variation de la distribution des diamètres. 
En particulier, il est bien sûr possible d'augmenter la taille des billes (au nombre de x N{) qui 
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ne portent aucune force. De telles approches discrètes permettraient d'explorer le sens précis de 
diverses notions (« classe granulaire dominante », « effet de paroi », « effet de desserrement »...) 
introduites dans le modèle de de Larrard [QQ|, utilisé pour la formulation des ciments et des 
bétons. 

5 Conclusion et perspectives 

Nous avons établi que la méthode de simulation utilisée permet d'obtenir de bons échan- 
tillons représentatifs pour les mélanges binaires considérés. La détermination de la compacité 
est relativement aisée, même pour des échantillons de quelques centaines de grains, lorsque 
le nombre de gros est de l'ordre de la dizaine. La compacité, pour le mélange binaire comme 
dans le cas monodispersé, augmente lorsque le processus d'assemblage comporte une phase 
d'agitation prolongée. Cet effet, bien que relativement modéré, se dégage nettement des marges 
d'incertitude statistiques. La compacité maximale usuelle correspond alors à la limite des as- 
semblages rapides. 

L'étude des état de compacité maximale peut être poursuivie et généralisée à diverses gra- 
nulométries, en particulier grâce aux performances croissantes des ordinateurs et au calcul pa- 
rallèle (pour lequel le LMSGC s'est doté d'une machine spécifique). Il sera particulièrement 
intéressant de relier certaines valeurs caractéristiques de la compacité (valeur maximale, valeur 
critique...) à diverses propriétés rhéologiques (dilatance, viscosité, frottement interne), éventuel- 
lement couplées à la ségrégation, des assemblages de grains seuls ou en présence d'un fluide. La 
simulation numérique des suspensions et des pâtes granulaires (suspensions très concentrées) 
est l'objet de recherches en cours au LMSGC (A. Lemaître, A. Hammouti, J.-N. Roux, P.-E. 
Peyneau) et au LCPC-Nantes (N. Roquet). 
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